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1. Introduction 

Matrix models describing the behaviour of Dp-branes originate in the observation of 
Witten that the massless modes propagating along the world volume of iV coincident 
D-branes are those of the supersymmetric Yang-Mills theory, obtained by the dimensional 
reductions of the d = 10 M = 1 theory down to p + 1 space-time dimensions. 

In various compactifications of string theory one encounters the nearly massless non- 
perturbative particles, obtained by wrapping the Dp-branes around vanishing p-cycles 
inside the internal Calabi-Yau manifold. Even in ten-dimensional Type IIA string theory 
there are solitonic particles , which are represented by certain black holes in the effective 
supergravity and are interpreted as Kaluza-Klein modes of the graviton multiplet in the 
compactification of M-theory on a circle Of course, these particles are no longer 

massless. 

Despite the variety of mechanisms by which such objects appear, their internal de- 
scription at low energies proves to be rather simple. In fact, if N such particles in d 
space-time dimensions are close to one another then their dynamics is described by the 
dimensional reduction of M = 1 super- Yang-Mills theory from d down to + 1 dimensions 
(first studied a long time ago for different reasons in j5||). The degrees of freedom in such 
quantum mechanics are represented by U(N) matrices X 1 , i = 1, . . . , d, together with the 
gauge field At and their fermionic partners. 

Although the exact computations in quantum mechanics of interacting particles are 
rarely possible, the supersymmetry allows one to get some exact answers. In this paper 
we shall concentrate on the correlation functions of the light-like coordinate operator. To 
state more precisely what we mean by that, let us consider the quantum mechanics with 
periodic time t ~ t+ 2n[3 and with periodic boundary conditions on fermions. In this case 
one can show that the observable 

R = Tt R P^f dt( A t+ X>) 

commutes with some of the supercharges (of course the choice of X 3 is arbitrary). In the 
limit f3 — > (and after Wick rotation) the computations in the quantum mechanics reduce 
to the finite-dimensional integrals, where A t becomes the 0-th matrix Xq = —iX^. Then 
the observables On can be expanded in 



Tr(X+)\ X+ = X 3 + iX\ 
2 



The paper is organized as follows. We are going to study the case d = 4 in great detail. 
We derive the determinant representation for the regularized generating function of the 
correlators of Tr(X + f and show that it obeys Hirota bilinear identities (when working 
with fixed chemical potential, e.g. in the grand canonical ensemble). Then we concentrate 
on the operators (Tr (X + ) 2 ^j and derive the asymptotics for the generating function in 
certain limits. We then briefly discuss d = 6, 10 cases. Then we proceed with the direct 
attack on the d = 4 integral for fixed but large N, using the saddle-point techniques, and 
derive interesting asymptotics both in the strong and in the weak coupling limit. In the 
weak coupling limit, we obtain agreement with the planar graph expansion. In the strong 
't Hooft coupling limit we get the agreement with the predictions from KP hierarchy. In 
the first appendix we check the strong coupling asymptotics by the direct quasiclassical 
calculation of the matrix integral. In the second appendix we demonstrate our method 
based on the KP differential equation on the example of the usual one matrix model. 

In the bulk of the paper we use the notation = X + , = X~ . We also denote by 
Z, F = logZ the partition function and the free energy at fixed particle number N and by 
Z, T = \ogZ the corresponding quantities at the fixed chemical potential fx. 



2. Supersymmetric matrix integrals 

2.1. Theory with four supercharges 

The dimensional reduction of the M = 1 SYM from d = 4 dimensions down to zero 
dimensions would produce a matrix model with 4 bosonic matrices X^, (i = 1,2,3,4 and 
2 complex fermionic matrices A , a = 1,2. All matrices are in the adjoint representation 
of the gauge group G, which we will take to be either U(N) or SU(N)/W. N . The matrix 
integral has the form: 

— ^ J dXdXexp ^ ^Tr[X M , X.] 2 +Tr\ & <r?[X lt , Ajj , (2.1) 

where cr M = (1, cr^), % = 1,2, 3, are the Pauli matrices. The two complex fermions A can be 
viewed as four real fermions, which we denote by x, rj, tp a , a = 1, 2: 

Ai = \ (77 - ix) » A 2 = \ + #2) 

and 

A„aa \ * 
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We also redefine the bosonic matrices as: 



-L(X 3 + zX 4 ), 0=-^=(X 3 - i X,) (2.2) 



and introduce an auxiliary bosonic field H (also in the adjoint). Then the integral in ( f2.1] ) 
becomes: 

/dX ce d^fj a dxdHd(J)d'nd(J) . . 
VoKG) 6XP( ^ } (2.3) 

S = (iTrHs + \TxH 2 + Tr[X Q , 0] [X a , 0] + §Tr[0, </>] 2 + . . .) 

where s = [X±, X2] and . . . represent the fermionic terms that are reconstructed using the 
following nilpotent symmetry of fl2.3|) : 

5X a = ip a , 5ip a = [(f), X a ] 
Sep = v , 8r)= [</>, 4>] 
S X = H : 6H=[<f>,x] 
5<p = 0. 

The symmetry 5 squares to the gauge transformation generated by (p, hence it is nilpotent 
on the gauge-invariant quantities. This symmetry was formally studied in in order to 
apply it to the model of ; it was powerfully exploited in |§ in the problem of computing 
the Witten index in certain quantum mechanical systems (first studied in the two-particle 
case in [0-^0], see also [llj). The action of the matrix integral (|2.3|) is 5-exact; in fact it 
may be written as 

S 



6 (-iTrxs-lTrxH-J^^MX*,® ~ |Trr#,0]j . (2.5) 



Now we proceed to reducing the integral (|2.3| ) to an integral with respect to the single 
matrix variable (p. The strategy is known for some time |12[ , and it consists of two steps. 

If the action is perturbed by the expression 8(. . .), which has a nice behaviour at 
infinity, then the integral should not change, which can be shown by doing an integration 
by parts. Consider the modification of the action S by the term 

S -> S + iSR, 

with 

R = ^e^Tr^Xp + k 2 ^x4>- (2.6) 



This perturbation makes the integral (2J.) localized near the zeros of H, <f>, x, tj in the limit 
of large k% : which can be shown by the saddle-point approximation. It reduces the integral 
( PTT| ) to a simpler one 

The behaviour of the integrand at large values of (p is still not good enough. To make 
it better behaved we modify the transformation 5. The current 5 is designed to respect 
the ordinary gauge invariance. In particular 5 2 is a gauge transformation generated by 4>. 
We wish to invoke yet another symmetry of the integral in ( |2.1| ) , which is the global group 
U(l) £ acting on the matrices X a , if; a via the rotations: 

e ie : X t + iX 2 i-> e ie {X x + iX 2 ) . (2.8) 

The rest of the fields are invariant under this U(l) £ group action. Let us denote the 
generator of this group by e. Then the new supercharge 5 acts as follows: 

5X a = </> Q , 8ijj a = [0, X a ] + ie ■ e^Xp 
5(p = r), 5r]= [0, cj>] 

r , (2 - 9) 

5 X = H, 6H = [(j>, X ] 
5<p = 0. 

The integral ( |2.1| ) has another U(l) symmetry (called the ghost number C7(l) g h) under 
which 5 has charge +1, the bosons X a , H have charge 0, the fermions %l) a have charge +1, 
the fermions x, r\ have charge —1, and the bosons <p an d 4> have charges +2 and —2, 
respectively. The measure and the action have the over-all charge 0. The modification 
( |2.9|) is consistent with the ghost number symmetry iff the generator e is assigned the 
ghost number 2. If we compute the modification of the action ( |2.7| ), we obtain a better 
behaved integral 

/ d %oh (G) a 6XP KlTT ^ X ^ X ^ ~ i< X " + + ^2) ■ (2.10) 

The factor Ki can be now reabsorbed into the X's and ifj's without affecting the measure, 
and the i/j's can then be integrated out. Also the matrices Xi,X 2 can be integrated out, 
producing the determinant 



Z(N - V) =1 Volfg) Bet < 2 - n » 
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The supersymmetry 5 allows a modification of the action by the observables 

S^S + Y^TnOr, 

n 

where SO n = 0. In our case the operators O n are simply the gauge-invariant functions 
of (p as they are also £7(l) e -invariant. The simplest operators whose correlation functions 
may be evaluated are the gauge-invariant functionals of </>, such as Tr(</>'). 

To summarize, we have shown that the computation of the (regularized) correlation 
functions of the observables Tr(</> n ) in the supersymmetric matrix integral (|2.1|) may be 
reduced to the computations of the integral over a single matrix <p of the form V(z) = 
-E n T n z n : 

-TWO) 

(2.12) 



Vol(G) Det (ad{(/>) + e) ' 
Remarks. 

1. In ref. ||, the similar perturbation has been used in the computations of the Witten 
index, which can be reduced to the computation of the integral ( |2.1| ) for the group 
SU(N)/1Z.n and without insertion of any observables. In that case the result of the 
computation was e-independent. Also, the integral over the eigenvalues of the matrix 
<fi was to be understood as a contour integral, to avoid the contribution of the flat 
directions that corresponded to the unbound free particles. In our case, the flat 
directions contribute to the correlation functions as well, and the parameter e serves 



as a regulator as in the computations of [13 



2. One may wonder about the physical meaning of the e-regularized integrals. Here it is: 

Z(N, e, V) = Tr w (-) F e- /3H e- eJ e- TrV ^, (2.13) 

where the trace is taken over the Hilbert space TC of the quantum mechanical system, 
H is the Hamiltonian, F is the fermion number, J is the generator of the global 
symmetry group (which we take to be SO(d — 2) for d = 4, 6 and SO (6) for d = 10, 



see below). For example, in the case d = 4, J = 2zTr(0[X L , X 2 ])- Just as in |1C 
this expression is related to the matrix integral in the f3 —> limit. One can also 
consider directly the quantum mechanical path integral, i.e. the integral over the 



space of loops. In this case the rational functions in the formulas ( |2.12|) and the 



similar formulas below are replaced by their trigonometric counterparts. Also one can 
consider a 1 + 1 model (Matrix strings) in which case the ratio of determinants leads 
to elliptic functions, just as in |Tj 
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2.2. Theory with eight supercharges 

This is the model obtained by the dimensional reduction of M = 1, d = 6 theory. In 
this model the index a of the matrices X a , if> a runs from 1 to 4. The symmetry U(l) e 
is extended to SO (4) « SU(2)l x SU(2)r. The matrices X a , form two copies of the 
representation (§, o) or * this group. Also, the fermion \ is promoted to a triplet x, which 
is in (i, 0). The same metamorphosis is experienced by the auxiliary field H — > H. The 
action is constructed by the same rules, the only difference being that 

where 

Si = [X^jXi] + \eijk[Xj, Xk\. 
The modification of the supercharge (|2.9|) is achieved by introduction of the generators e_L© 



tR = ( £l 2^ £2 ) © ( £l 2 62 ) °^ the Cartan subalgebra of S'0(4). The modified transformations 



are: 



Si(ji = [(f), X{\ + ieiX 2 , Stf>3 = [(f>, X 3 ] + ie 2 X 4: 

5ip 2 = [<P, X 2 ] - ie]Xi, £0 4 = [0, Xt] - ie 2 X 3 

5 X i = Hi 5H s = [(f>, X 3} 

8Hi = [(f), xi] + 2itLX2 SH 2 = [(f>, X2] - 2ie L Xi- 



(2.14) 



Now we get, instead of ( 2.12|) , the following one-matrix integral 

f d(f)e-^ v ^ Det (ad((f)) + ei + e 2 ) 

J Vol(G) Det(ad(0) + ei)Det(ad(0) + e 2 )' ' 
£.5. Theory with sixteen supercharges 

It is of great interest to obtain the similar expression for the integrals occurring in 
the reductions of d = 10 SYM. In this theory the matrices X a ,if> a have the index a 
transforming by the 8 of the group SO(8). The antighost x belongs to 1 © 6 of SU(4) C 
SO(8). Introduce the notation 

B, = X 2l _ x + iX 2l ,i = 1,2,3,4. 

The matrices Bi are in 4 of SU(4) and B\ are in 4. The "gauge condition" s splits as: 

4 

(2.16) 



$y = [fl i ,.B J -]+Jey W [fl fc ,.B I ], 
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The action constructed by the standard rules coincides with that of the dimensional 
reduction of d = 10, M = 1 SYM. The gauge field has ten components, which become 
0, (p and X a . The sixteen-component fermion splits as ijj a , with eight components, x with 
seven components and r\. 

The global group SU(4) (which is not to be confused with the .R-symmetry group 
of Af = 4 SYM in four dimensions!) as before allows us to modify the supercharge 8 in 
a manner analogous to (|2.9| )-( |2~T4"| ). The Cartan ganarator of SU(4) may be written as 
ei © e 2 © e 3 © (e 4 = -e x - e 2 - e 3 ). The integrals ( |2TT2l )-( |2TT5| ) generalize to: 

j. p -TrVW 



Vol < G > (2.17) 
Pet (ad((p) + ei + eg) Pet (od(0) + e 2 + £3) Pet (ad((p) + e 3 + gj 

Pet (ad(4>) + ei) Pet (ad(0) + e 2 ) Pet (ad(0) + e 3 ) Pet (ad(<f>) + e 4 ) ' 



3. Determinant representation of the correlation functions in the d = 4 case 

In this section we study in detail the grand partition function 



Z(p,e,V)=Y i e* N Z(N t e,V). 



N=0 



We show that Z(fi, e, V) has a determinant representation, very much like the correla- 
tion functions in Sine-Gordon, and related models are expressed in terms of Fredholm 
determinants 111511 . 



3.1. Eigenvalue integral 

First we write the integral ( |2.12| ) in terms of the eigenvalues icpi, z^jv of the anti- 
Hermitian matrix cf>: 

Z(N,e,V)=f ^^n^i^ne-n*), (3-1) 
J u n A^!(27re) iV (pi - fa + it - L .- L 



where we changed V(ix) to V(x). The integral ( p.l|) can be rewritten, using the Cauchy 
formula, as 



It turns out that the grand partition function (that is, with fixed chemical potential) 
can be written as a Fredholm determinant of an integral operator. Let us introduce the 
notation 

w '^ v 1= /..ntr .-l.-....v *w s *>- («) 



A = A 2tt e-i(x k - Xk+i) 



We may rewrite the sum over all the elements of the permutation group in (|3.2| ) as the 
sum over the conjugacy classes that are labelled by the partitions of iV: 

oo 

N = J2 l dh di>0. 
1=1 

Every permutation in the conjugacy class, labelled by d = (d±, d%, . . .), is similar to the 
product of cycles of lengths 1,2,..., the number of times the cycles with length / appear 
being precisely di . The number of permutations in the given conjugacy class d is equal to 

N\ 



W dl d, 



and the sign of any permutation in this class is (— ) ! (— ) N . Thus, ( |3.2[ ) may be repre- 
sented as 



(3.4) 



and the grand partition function is equal to 



^, e ,y)=exp^(-y-V-^l^. (3.5) 
I 

The quantity Wi(e, V) may be represented as Wi = TtK 1 , where K is a linear operator 
acting in the space of functions of one variable, as follows: 

(*/)(*) = e-™ / £ f } V) y (3.6) 
Jr 2tt e - i(x - y) 

Therefore, the grand partition function becomes 

f_y-i 

Z(fjL, e, V) = exp —, — Tr(e M -ftT)' = Det (/ + e"K) . (3.7) 
I 



3.2. Another representation for quadratic V 

Let us consider the case of a Gaussian potential V(x) = | (jTj an< ^ wr ite the parti- 



tion function ( |3.1] ) again as a matrix integral 

Z(N, = J ^fy exp (-iTr[0, X] 2 + ^Tr</> 2 - I e 2 TrX 2 ) . (3.8) 

Considering the matrices X and <f> as the Hermitian and anti-Hermitian part of the same 
complex matrix Z = (e/y/QX + (\/£/e)</>, we rewrite the matrix integral as 

Z(N, = J dZalZ^ e" s , S = §Tr[Z, Z^f + ^TvZZl (3.9) 

In polar coordinates 

Z = UH 1 / 2 

where U is unitary and if is a Hermitian matrix with positive eigenvalues y±, . . . , y^, the 
measure and the action read 

dZdZ ] = dUdH, S = TrH 2 - TrU^HUH + ±£TrH. 



Using the Harish-Chandra-Itzykson-Zuber formula PB][I7| we perform the ^-integration 
and find 



1 ' ]N 



Z ( N >® = m'-- Ve 'Dety(e™")- ( 3 - 10 ) 



The grand partition function will be expressed in terms of the quantities 



as 



r 1 

Wt = Udyi e -Mfw+(»«-»H-i) a ] j ( w = (3.ii) 
Z(/U) = exp V<^. (3.12) 



i=i 



Clearly = Tr/C', where K, now acts on the functions on the positive semi-axis as 



/ 1 

1Cf{y) = e~^ y / e-^v-y') f{ y ')dy'. (3.13) 
Jo 

Therefore we arrive at the same determinant representation ( |3.7| ) where the kernel K is 
replaced by its Fourier transform /C. 
For small 



'-y-'Wt = ± gi (0, (3.14) 



where gi is an analytic function. 
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4. The grand partition function as a tau-function 

In this section we show that the grand partition function 



~ r 00 N n (x -x ) 2 

can be represented as a tau-function of the KP hierarchy. 

4-1. Vertex operator construction - bosonic representation 
Introduce the bosonic field <p(z) with mode expansion 

[Jn, Jm\ = n$m+n,0; [£, <?] = 1 (4-3) 

and the vacuum state |/) defined by 

J n \l) = 0, (n>0); p|Z) = Z|0- (4-4) 

The associated normal ordering is defined by putting J n with n > to the right. Define 
the vertex operator 

V e (z) =: e ^ (2+ie/2) :: e ~<p(z-ie/2) : , (4.5) 

which satisfies the OPE 

Ve(z)V e (z') = {z ^~pl e2 ■ Ve(z)V e (z') : , (4.6) 

the Hamiltonian 

#[£] = ^£ n J n , (4.7) 

n>0 

and the operator 

^ = exp(e^| ° ^V e (z)). (4.8) 
Then the vacuum expectation value 

r [t) = (0|e H t%|0> (4.9) 
11 



is equal to the canonical partition function ( f4.1|) , with chemical potential \i and potential 
V(z) =u(z+ l -e)-u(z- l -e)=-jt 



2 J \ 

oo 

U(z) = -J2t 



oo 

n 



n=0 (4.10) 



n=l 



where can be reabsorbed into the definition of V and where the coefficients T n and t n 
are related by 

r «-i = ^H fc KVw2, (4.11) 



4& 

fc=0 



Fermionic representation 



The fermionic representation of the partition function is constructed using the 
bosonization formulas 

^(z) =: e-vW : ij)*(z) =: e^ z) : =: V*(*)V>(*) : , ( 4 - 12 ) 
where the fermion operators 

V(z)= E Vr«- p -i, V*W= E ^ _r_i (4-13) 

satisfy the anticommutation relations 

[^,V:]+ = 5 rs . (4.14) 
The operators (|4.7|) and (|4.8| ) are represented by 

H[t} = J2 t nY, : W-n^r-- (n € Z) (4.15) 

n>0 r 

r r°° dx e el 

At =expe /x / — : ^(a; + - i-) : . (4.16) 

L J-oo 2m 2 2 J 

and the vacuum states with given electric charge / satisfy 



(l\^_ r = <Z|W* = (r>Z) 
^10=^10 = (r>0- 

12 



(4.17) 



The original expression (O) is obtained from the expectation value (|4.9| ) by first commut- 
ing the operator e H to the right until it hits the right vacuum by using the formulas 



(4.18) 



and then applying the Wick theorem to calculate the expectation value 

ffln^wio=nfe)'n j;:ga,~-g - w 



4-3. The KP hierarchy 



The partition function (4.9) is a particular case of the "general solution" of the KP 



hierarchy obtained as the limit N — > oo of a general iV-soliton solution [|18 

n[t} = (l\e^>o tnJ -n a \l), (4.20) 

where the GL(oo) rotation 

fi a = exp ( J dxdy a(x,y)ip(x)ip*(y)j (4-21) 

is parametrized by an arbitrary integrable function a(x, y)0. 
The tau- functions 77, 1g2, are Fredholm determinants 

r ; = det(l + e^Ki) (4.22) 

of the kernels ^ 

jjf,(a;,y) = If + U ] E t {x) = x l exp ( Vt n / ) . (4.23) 

— ze) x — y — lie \ ) 



For / = we get precisely the operator K (|3l]). 

The KP hierarchy of differential equations is generated by the Hirota bilinear equations 



y^^'exp(^(t n -0^)r^t n -V^ r v {t' n + h-^=Q (/'</). (4.24) 



n>0 



1 In the soliton solutions it is a sum of delta- functions, a(x,y) = y\ ai5(x — pi)5(y — qi), so 
that the operator f2 a is a finite sum, fi a = y\ —1 aii^(pi)ip* (qi). 
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Let us sketch the proof of ( 4.24 ). First we remark that each element G GL(oo) is 
represented by an infinite c-number matrix a = {a rs } 



o^fi -1 = ^^ s a sr , n- 1 tp;n = J2 a rs^t- (4.25) 



As a consequence, there exists a tensor Casimir operator 

dz 



.V(*)®V*(*) ( 4 - 2 6) 
2m 



r 

which satisfies, for any Q G GL(oo), 

S12 O = n Si 2 . (4.27) 

On the other hand £12 |Z) <8> |Z) = because, according to ( f4.17| ), for each r either ip r or ?/;* 
is annihilated by the right vacuum |/). Therefore ( |4.27| ) implies that Si 2 0|Z) (g> f2|Z) = 0. 
Taking the scalar product with (I + l\e H ^ ® (Z — l|e H ^ ' we find 

f^Ti {l + ^M^WM " = 0, (4.28) 

where the integration contour surrounds the origin. Equation Q4.28| ) simply reflects the 
fact that the tensor Casimir ( [4.26| ) is constant on the orbits of GL(oo). Finally we use the 
bosonization formulas ( |4.12| ) to represent the fermions as vertex operators, i(j(z) — > V_(z), 



i(j*(z) — > V+(,z), where V±(#) act in the space of the coupling constants as 



V ± (,) = exp ( ± £ t „z») exp ( T In T £ £_ ) . (4.29) 

n=0 n=l 

The general case Z 7^ I' is treated similarly, and one obtains the following identity 

jdz (v+{z) • 77[i]) (v_(s) • T V [t']) =0 (/' < Z), (4.30) 

which is identical to the Hirota equation ( f4.24| ) . 

The differential equations of the KP hierarchy are obtained by expanding (|4.24|) in 
the differences y n = \{t n —t' n ). In the case V = I, the first non-trivial equation (the KP 
equation) is obtained by requiring that yf term vanishes: 

(£t + i- 4 iiB T,[t+y]T,[t - y] l--r - (4 - 31) 
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In terms of the "specific heat" 

Q2 



dt\ 

the KP equation reads 



«M = 2 ^2 lo g^ (4-32) 



d 2 u d 



du „ du d u 



4^ h 6u- h 



<9t 3 9ti 9tf 



0. (4.33) 



5. The d = 4 integral with quadratic potential: KP equation, weak coupling 
and double scaling 

In this section we study the case t n = 0, n > 3. 

5. i. Reduction to a single equation 

A potential of the form V(x) = —fi + Ax 2 , A = £ 2 is related to the three couplings 
ti,t 2 ,h by ^ 

= ie - ~e 2 £ 3 - |0 , £ 2 = -3kt 3 . (5.1) 

Rescale e — > 1. Then w = —2d 2 J 7 and it is easy to show that Q4.33|) implies the following 
partial differential equation for the function %j) = £it(/i, £): 

^ + V^m + | - = (5-2) 

By comparing it with the expansion in (|3.12| )-( |3.14j ), we see that a(£) = 0. If we expand 

00 2 

^(^0 = ^e^e-^(' 3 -S(0, (5-3) 



then eq.( |5.2| ) is equivalent to the infinite system of recursive first-order differential equa- 
tions: 

vfi= l 2 E W^ Pq - ( 5 - 4 ) 

p+q=l 

The form ( |5.3|) is dictated by the semi-classical approximation to the integral ( |2.10|) . 
Indeed, the expression j^(l 3 — l) is nothing but the classical action evaluated on the solution 
to the equations of motion: 

[X 1 ,X 2 ] = -2i\<f>, 
[ ( />,X 2 } = -iX 1 , (5.5) 
[<j ) ,X 1 ] = +iX 2 . 
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The solutions to (|5.5| ) are classified by the decompositions of iV-dimensional representation 
into irreducibles of 577(2). The logarithm of the grand partition function takes into account 
only irreducible /-dimensional representations, the rest is generated by the exponentiation. 
The functions rji therefore describe the quantum fluctuations around the saddle points. 

We conclude this section by listing the two equivalent forms of the equation obeyed 
by u: 

fa + W>n + |0 - VwOm = 

1 1 (5 - 6) 

2u\ + + uu^ + -(u - u^)^ = 0. 

5.2. Weak coupling limit 

For low values of Z, eqs. ( |5.4| ) can be solved explicitly. It is interesting to look at the 
large A = £ 2 asymptotics of the solutions. We expect that as A — > oo the partition function 
at fixed N has the following scaling behaviour of Z(N, A) = e F ( N >^>: 

Z ( N > A ) = (1 + -..) = (1 + •••) ■ ( 5 - 7 ) 

On the other hand it follows from fl5.4f ) that 

6 m 6 2m- (2n-l)!! 

^ = -y + 7^E c2n + ••• ( 5 - 8 ) 

and therefore indeed 

Z(1,A) = 1, 2(2,A) = -lg=^, (5.9) 

in accordance with the scaling ( |5.7| ). 

It turns out that equation (|5.2|) has another interesting property. Suppose we are 
studying the 't Hooft limit, where the free energy has an expansion of the form 

oo 

F(N, A) = N 2 - 29 F g (N/X). (5.10) 

Given Z{ji, A) = e^^^' we extract the fixed- N partition function via the Fourier transform 

Z(N, A) = j ^L e -^+Hi^) , (5.11) 
16 



which can be taken, in the large-iV limit, using the saddle-point approximation. In the 
planar limit (N = oo, with X/N finite) the functions A) and F(N, A) are Legendre 
transforms of each other. If we keep all the 1/N corrections, then the corresponding 
expansion of A) is: 

oo 

(5.12) 



and 



u 



5=0 



(/*, A) = u 9 (MA) A" 25 , u 9 (x) = -2J*(x). 
g=o 



Let us introduce the variables 



(5.13) 



x = /i/X, y = X 2 , x(x,y) = u(fjL, A), 



(5.14) 



which are relevant for the 't Hooft limit. Then eq. (|5.2| ) may be rewritten as: 
X - 2 fx - ^ J - 4yxy + XXx - = 0, 



(5.15) 



which after the expansion 



9=0 



X) 



reduces to the infinite system of recursive equations for x 3 's 
1 



--2x + xo)x' g + (±-4g-Xo) = \ 1 - 

a-\-b=g ,a,b^Q 

Xo ~ 2 - J xo + XoXo = 0- 



The equation for xo is the only non-linear one. Its solution is: 



XaX b , 9 > 



(5.16) 



(5.17) 



«0 2 1 

x = Y Xo + Xo + T2- 



(5.18) 



Of course eq. (|5.2j ) has more general solutions, in particular those for which the 
expansion ( |5.16| ) is not bounded as g > 0. It turns out that the solution corresponding to 



the matrix integral in question does have the form (|5.16| ). In Appendix A we show that 

«0 = — 7T 2 . 
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It follows that, for large iV and A, the free energy has 't Hooft-like behaviourl: 

5.3. Double scaling limit near the quadratic singularity in x( x ) 

So far we investigated ( |3.1|) only in the large- ('t Hooft) limit in the canonical 
ensemble or, equivalently, in the large-/U limit for the grand canonical ensemble. 

As we know from and Appendix B, the universal scaling behaviour of higher 1/N 
corrections sometimes can be summed up to some functions obeying non-linear differential 
equations, such as Painleve I for the pure 2d gravity. 

It is reasonable to ask whether we can do the same with the 1/// expansion for our 
model starting from the general KP equation ( |5.15|) and what the physical or geometrical 
meaning of this expansion is (we recall that in the pure gravity described by the one-matrix 
model of Appendix B the corresponding 1/N expansion has the meaning of the expansion 
over the genera of the topologies of the two-dimensional manifold) . 

Let us concentrate on the square-root singularity of ( |5.18| ) at Xc = ^> x c = — + 
We try the following ansatz: 

X = Xc + y a v(z), z = y b (x-x c ). (5.20) 

As in the case of the one-matrix model, the presence of a quadratic singularity implies 
that b = —2a. In full analogy with Appendix B (the only difference being that y takes the 
place of 1/N 2 ), inserting this ansatz into (|5.15|) , and neglecting the sub leading corrections, 
we obtain b = — I, o = | and that function v(z) satisfies the Painleve I equation: 
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v" - 3v 2 + -=z = 0. (5.21) 

From this equation we find the following coefficients of the 1/A expansion in ( |5.12j ) 
(which is the same as the 1/fj, expansion in this approximation) for the singular part of 
the free energy near the critical point: 

J- sing = -At, ^ ing = -^logA, (5.22) 



2 Recall that F = T — uN in our conventions. 
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where A = const (x — x c ) (we choose the constant in such a way that the coefficient in front 
of A 5 / 2 be —1, then the next coefficients, — • • •, are universal constants). 

So everything goes just as in the pure 2d quantum gravity. The l/fi expansion looks 
like the topological 1/N expansion, the coefficients giving the leading scaling behaviour of 
the partition functions of successive topologies (see the details in [|0|)- ^ i s tempting to 
speculate that the quadratic singularity in x( x ) corresponds to the pure gravity. It may 
be related to the large planar graph expansion with respect to g in the model of dense 
self-avoiding random paths (we will argue at the end of section 6 that our matrix integral 
describes such a model in the large- N limit). It would be interesting to demonstrate it by 
passing from the grand canonical to the canonical ensemble for the free energy. 



6. Saddle-point approach 

6.1. d = 4 integral 

So far we managed to calculate the grand canonical version of the integral ( |3.5|) in the 
large-// limit by the use of the KP equations. It is not clear whether we can derive from this 
asymptotics the large N limit of canonical partition function. In fact, we shall show that 
it is possible by comparing the results with a more direct approach, originally proposed in 
this context by Hoppe Namely, in the case of Gaussian potential V(x) = Xx 2 , it is 
possible to solve the integral saddle-point equation for the distribution of the eigenvalues of 
a matrix in ( |3.1| ) . We work out the details of the solution (correcting some minor mistakes 
in [2T and actually deriving the result) and extract interesting critical behaviours of our 
system. 

It is natural to scale the coupling A as N, and rescale e to 1, i.e. to set: 

A=-, e-1. 6.1 
9 2 

Indeed by rewriting the integral (|3.1| ) as: 

Z(N,X) = 1 —^ T [ dXdcf>exp f-Tr (</> 2 + X 2 ) + ?-Tr[X, cj)] 2 ) (6.2) 

vol(G)A— J V A / 

we see that A= plays the role of the coupling constant while y is the 't Hooft coupling. In 
the large- N limit the integral ( |3.1| ) localizes onto the critical point of the effective potential: 

VM) = V^i) + E log (l + (0 ._^. )2 ) • ( 6 - 3 ) 
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Its critical point is found from the equation: 



1 



N W0 fc -0 .)(i + (0 fc -0 .)2)- 

3r K 



(6.4) 



In the usual fashion we assume that the eigenvalues in the large- N limit form a continuous 
medium of density 

P( x ) = ^H 6 ( x - ( t>i)- ( 6 - 5 ) 

i 

When g is real, it is natural to expect p to vanish outside of the interval [—a, +a], and to 
be an even function p(x) = p(—x). We introduce the resolvent: 



W{z) 



z — x 



(6.6) 



and rewrite the eq. (|6.4|) as: 



2x 



— = W(x + iO) + W(x - iO) - W(x + W(x - i) 



(6.7) 



and where it is assumed that x G [—a, a]. The spectral density is given by the discontinuity 
of W(z) along the cut at (—a, +a): 



W(x + iO) - W(x - iO) = -2yrzp(x). 



(6.8) 



Let us introduce the holomorphic function 

~/2 



G(z) = -x + i 
9 



(6.9) 



which is related to the expectation value of the fermionic current (^*(z)i(j(z)) in the grand 
canonical ensemble ( |4.1| ). The saddle-point equation (|6.7| ) can now be rewritten as 



G [ x + 



G x 



The definitions 



and ( |6.9| ) imply that, in case of positive coupling A = N / g 2 



(6.10) 



W(z) = W(z), W(z) = -W(-z), 



G(z) = G(z), G(z) = G(-z) 



(6.11) 



and also that the function G(z) has the cuts at (±| — o,±| + a). It is also clear from 
( |6.10|) and (|6.11| ) that G(z) is real when z G IR, iM, (±| — a, ±| + a). Hence G(z) defines 
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a holomorphic map of the region U bounded by IR_|_, zIR_|_ and by the sides of the interval 
(|, | + a) onto the upper half- plane 7i. The inverse map is given by the following integral 
formula: G(z) = £, where: 



z = A 



dt(t — xs) 



51 ~ Xi)(t ~ X 2 )(t - £4) 

The map acts on the special points X\ > X2 > x 3 > x 4 and 00 as follows: 



c 


z 


+00 


+00 


X\ 





x 2 


i 
2 


x 3 


| + a 


x 4 


i 
2 


— 00 


-Moo 



These conditions imply the following equations on Xi, A, a: 

t-X3 



\=A dt 

J x 2 

a = A I dt 



' x 3 

a = A dt 
J X4 



V(t- 


x 2 ){x 1 


-t)(t 


- x 4 ) 




x 3 - 


- 1 






- t)(xi 


-t)(t 


- x 4 ) 




X3 - 


- 1 






- t){x\ 


-t)(t 


- Xa) 



From ( |6.9| ) we know the large- z asymptotics of £: 

1 



c 



;Z Z +Z~ Z +5Z-* + .. 



where 



A 1 A _3 5 5 



1 /"^ a 

S = — + v ~ ^..,.,2, 
4 



(6.12) 



(6.13) 



(6.14) 



p(y)y d v- 

At large t the function z(Q, as given by (|6.12|) , has the following expansion: 



(6.15) 



z = 2A I C 2 +a + ai( * + a 2 ( ' + a 3 ( i + 



(6.16) 
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where 

x 1 (t-x 3 ) 1 



y/(t - Xi)(t - x 2 )(t - x 4 ) Vt - XI 



ao = dt 

J X\ 
(~) k ~ 1 

ak 2k-l ^ k+x ^ k ~ 1 ^ k>0 (6 ' 17) 

7,= E ("/) (7) 

By comparing ( |6.14| ) and ( |6.17| ), we get the following equations on Xi, A: 

A = 9 - 
2 



a (xi) = 

X\ + X2 + X4 = 2X3 
2|2|2_r)2 _ ^ 



(6.18) 



which, together with eqs. ( |6.13| ), fixes everything completely. As shown in [21], the 
equation ao = follows from ( |6.13|) by a contour-deformation argument. We introduce 
more notation: 



, V% 
Ui = gxi, Xi 



yi - 2/4 

y^-yi 1 ^ ^ n '1 1 

m = , 1 >m > 0, m = 1 . 

2/1-2/4 m ' 

with which eqs. (|6.13| ) assume the following form: 

(2/2 + 2/4 - yi)K(m)+2(yi - 2/ 4 )E(m) = 



(6.19) 



-(2/1 + 2/2 - 2/4)K(m')+2(z, 2 - y 4 )E(m') = yj ^j^. 
where we use the standard elliptic functions 



(6.20) 



7T 7T 

K(m) = / =, E(m) = / d#Vl - msin 2 #, (6.21) 

./o VI - msin 2 6> ^0 



which have the following crucial properties: 



K(m') = v / mK(l - m), E(to') = -=E(1 - to) 

'TO 



E(m)K(l - to) + E(l - to)K(to) - K(to)K(1 - m) = -. 
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(6.22) 



In the sequel we use the short-hand notations E = E(m), K = K(m), d = E/K. The first 
equation in ( |6.20| ) allows us to express A2 in terms of m, while the second, together with 
( |6.22|) , gives g{y x - y 4 ): 



1 



A 2 = 1 - 20, g{ yi - 7/4) = — K . (6.23) 

From (|6.18| ) and the equations A4 = A2 — m and Ai = A4 + 1 we get an expression for 
yi — 2/4 and consequently for g: 



K 4 

9 2 M = — (4mA 2 + 1-3X 2 2 - 2A 2 ) 
K 4 

(-3tf 2 + 2(2 - ra)tf - (1 - m)) . 



(6.24) 



3tt 4 

We can now compute -^(Tr<^> 2 ) = f and F(N, g): 
9 4 dF(N,g) 



v{m) 



N 2 dg 2 



l_ / 4mA 2 (l - m) + (5A| - l)(2m - 1 - A 2 ) \ 

12 + 5tt 2 V 4mA 2 + 1 - 3A 2 - 2A 2 ^ 1 } 

1 K 2 10t? 2 (t9 + m — 2) + 2-&(6 — 6m + m 2 ) + (1 — m)(m — 2) 
~ 12 ~ 5^2 3t? 2 + 2(m -2)d + l -m ' 

The formulae ( |6.24j ) and (|6.25|) give the exact analytic solution of the large-iV model in 
the parametric form. 

Small-g expansion. In this case we expect to get a regular planar graph expansion 
of the matrix integral ( |3.8| ) with respect to the quartic term in the action. The careful 
analysis shows that the g — > limit corresponds to m — > 0. In this limit we can expand: 



2 1 2 1 3 119 4 

q = m H m H m + . . . 

y 128 128 16384 

1 o 1 q 133 a 
-m H m H -=m + . . . 



(6.26) 



which implies 



256 256 4096^ 



1 2 1 4 , 
V= 2 9 ~2 9 + 



F(N,g) =N 2 Qlog^ 2 -^ 2 + .. 
in perfect agreement with the planar graph expansion and formula ( p.7|) 
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(6.27) 



Large-g limit. This limit corresponds to the situation where we are far beyond the 
convergence radius of the ^-series. The integral (|3.8|) is dominated by the commutator 
term in the action and hence the fluctuations of the matrices are very large due to the 
presence of the zero modes. It follows from the careful study of (|6.23|) and ( |6.24|) that the 
g — > oo limit corresponds to m — > 1. In this limit (m = 1 — e): 



1 



9 



V 6 J 



12^ ' " > 



1 



/16 



v 



+ 2o^ log It 



iW(t)- 



(6.28) 



Hence 



(12tt)I 4 

- ff 3 



20 

F(N,g) = -N 



(12tt) 



3 a 
-# 3 



2 , 3(12tt)3 _2 



£ 3 +■■■ 



(6.29) 



40 5(12tt)1 

in perfect agreement with (|5.19| ) ill. In fact, the strong coupling expansion can be greatly 



simplified if we choose L 



log ¥ 



as an expansion parameter and systematically neglect all 



non-perturbative (in L) corrections, i.e. we consider the leading logarithmic approximation. 
We then get the very simple formulas : 



1 



9 



(1 - 3L) 



12tt 4 L 3 

1 /1-10L + 20L 2 \ I 



v = 



20n 2 L 2 



1-3L 



+ 



12' 



(6.30) 



The surprise is that formulas (|6.30|) are exactly equivalent to ( |5.18|) . 

2 

Proof. Equation (|5.18| ) with a = \ gives: 



T> - n v 3 1 v 2 T - 71 v 5 ^ v 4 + 1 V 3 

The Legendre transform leading from JF to -Fq yields: v = x — 2^, while # 2 
Clearly, it leads to (|6.30| ) if we substitute 



Xo 



7T 



(6.31) 



(6.32) 



3 Indeed, & 



- X :i _ 3(127r)5 
~ 40 
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This confirms once again that indeed ao = \. 
Remarks. 

1. It is very tempting to speculate that the relation to supersymmetric gauge theories, 
which was one of the original motivations of this work, is somehow revealed by the 
appearance of a family of elliptic curves, parametrized by the value of coupling g 



just as in ||22|| . Notice that the solution that we studied here has something to do 
with mirror symmetry Indeed, the naive coordinates in the space of our Lagrangians 
(which is just the coupling g for a quadratic potential) have been replaced by the 
period of a certain differential on the elliptic curve. It is conceivable that a similar 
construction takes place for more general potentials. 
2. The partition function does not change, in the large-iV limit, if we substitute the 
commutator by the anticommutator in the action in (|3.8| ). The latter model generates 
the statistical ensemble of </> 4 -type random graphs covered by dense non-oriented self- 
avoiding random loops and has been studied in P3|. It describes the dense phase of 
the 0(n) loop-gas model ||24j| with n = 1. It also can be viewed as a matrix model 
counting 3-colored planar graphs of the 4> 3 type: each propagator has one of 3 colors 
and all three propagators meeting at any vertex have different colors. This 3-coloring 
problem can be formulated as the following matrix model of three hermitean matrices 
B(lue), W(hite) and R(ed) with the action: 

S = NTr(gBWR + gBRW + B 2 + W 2 + R 2 ) (6.33) 

To see that it coincides with the original Hoppe's model it is enough to change the 
sign of one of 2 cubic terms (which is equivalent to the change of anticommutator of 
W, R to the commutator) and integrate over one of the matrices, say, B. The critical 
behaviour (thermo dynamical limit) is due to the dominance of graphs of infinite size, 
which renders the (/-expansion of the partition function divergent and is therefore 
determined by the singularity in g 2 closest to the origin. The latter is a solution of 
the equation g'(m) = 0. It should appear for negative g 2 and corresponds to the 
situation where all three cuts are located on the real axis, symmetrically with respect 
to the origin. When g increases, the end-points of the cuts get closer and a singularity 
occurs when they touch one another. One can show, using the symmetry ( |6.11| ) of 
W(z), that when g 2 is real and negative, the saddle-point equation (|6.7| ) becomes 
identical to eq. (3.3) of ref. fl2~5 |. It is known [2JJ that the critical behaviour of the 
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dense 0(1) model is in the universality class of the pure 2d quantum gravity. For 
example, the one-point function behaves as 

Here is the explicit formula for (g 2 )': 

= -^K^-^-\ +m \ (6.34) 
dm m(l — m) 

3. The last assertion can be partially confirmed by the study of the specific heat at 

the fixed chemical potential. From the eq. (|5.18 ) we get the critical point x c = 



— (^2 + ^2), which can be inserted into eqs. ( |6.30|) - (|6.32|) to yield g 2 = —3^4, which 
is negative indeed. Note, however, that the corresponding value of e = 2.16 is way 
larger than the leading log approximation allows us to see. 

7. On the correlation functions in the d = 6, 10 cases and directions for the 
future 

This section is devoted to a work in progress. We sketch the possible similar saddle- 
point approach to the d = 6 integral. We also attempt a fermionic representation for the 
d = 10 integral. 

7.1. Saddle-point approach to d = 6 integral 

We keep the same notations for the resolvent F and density p. We set e\ + €2 = 1, 
ei = /3, 62 = 7. Equation ( |6.7|) is replaced by: 

^ = F(x) - \ (F(x + i/3) + F(x - i/3)) 

9 (7.1) 

- \ (F(x + 17) + F(x - i 7 )) + \ {F{x + i) + F(x - i)) . 

Let us introduce the derived functions: 

m = ^- + F(x+ i -§.)-F(x- i -§.) 

U il ^7 X (7 ' 2) 

The saddle-point equation (|7.1| ) is equivalent to 

9 { X+ 2 i ) +9 { X ~2 i ) =0 " (7 - 3) 

The function g(x) has four cuts: at x G ±|, ±|(/3 — 7) + (—a, +a), it is real: g(z) — g(z) 
and it is purely imaginary for z G ilR. It would be nice to guess the correct function from 
the stated properties. We plan to return to this problem in the future. 
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7.2. Fermionic representation for the d = 10 integral 

We now proceed with the fermionic representation of the d = 10 integral ( |2.17|) . 
Unfortunately we were not able to find such a representation for all values of €1,62,63. 
However, let us consider the limit €3 — > — 6\, £4 — > —62- At the same time we keep 

e M (ei+e 3 ) =e^ 

finite. We claim that, in this limit, 

Z^ £ = (0\e H ^\0), (7.4) 

where 

r+°° dz 

Q u = e» : ib(z-aWz + aU*(z-bU*(z + b) : (7.5) 

J-00 2m 

with a = h(6i + £2)7 & = |(ei — ^2)- Now the relation between V and J7 is modified into: 

V(z) = U(z + a) + U(z - a) - U(z + b)- U(z - b). (7.6) 



As in the index computation it is possible that, by making the appropriate mass 
perturbation, we reduce the d = 10 integral for the gauge group U(N) to the products 
of d — 4 integrals for gauge groups U(ni) x . . . x U(rik) with N = YluLi ^ n i- Under this 
assumption: 

00 00 
Z d = 10 (^ 6, V) = n Z d= \l^ e, V) = J] Det( J + J»K) = Qfa K). (7.7) 

1=1 1=1 

The latter expression is very interesting, since it possesses certain modular properties 

and allows one to deduce the large- N asymptotics using very little information about the 
operator K itself: 

G(p, K) ~ exp 2^/N\U 2 (-K)l (7.8) 

where 

00 1 

Li 2 (-K)=^-Tr(-i^. (7.9) 

l=x 
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8. Conclusions 

Here we summarize the results of our computations. The integral ( |3.8| ), which is a 
cousin of ( |3.1|) , is studied in two regimes: at fixed N and at fixed fx. In the first case we 
got the large- N asymptotics in the 't Hooft limit, see (|6.24[) and ( |6.25| ). In the second 
case we showed that the grand partition function is a particular tau-function of the KP 
hierarchy. In particular, we obtained eq. (|5.2|) for the specific heat u = — 29^JF in the case 
of quadratic potential V ~ Xz 2 : 



1 1 

2u x + -u + uu^ + ~(u - u^jn = 0. 
A o 



In the large- fi, X limit, we obtained the simple explicit formula for the specific heat as a 
function of x = fx/X: 

-— u 2 + u + — = x. (8.1) 

We observed various similarities to the properties of supersymmetric gauge theories in four 
dimensions and we hope that our results will find their place in the study of dynamics of 
.D-particles in various dimensions. 
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Appendix A. Determination of cto 

We now present a trick that allows us to find the value of the unknown coefficient ao 
in ( |5.18|) . If ao 7^ then for large x one has u ~ ±w — and 



(A.1) 
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Below, we rescale A — > A/2 to be in agreement with the notations of (|3.13|) . 

Let e~ Ek be the eigenvalues of the integral operator /C. From the determinant repre- 
sentation (3/7) of the partition function, it follows that e^^'^ = Yl k (l + e fl ~ Ek ) . Corre- 
spondingly, the mean value of the number of particles is given by 

w = v=£ttpw- (a ' 2) 

k 

We are interested in the limit where both \i and £ (and therefore Ek, see below) are very 
large, i.e. a kind of low-temperature limit for a Fermi-gas with energy levels given by the 
spectrum of the operator log /C. In the low-temperature limit we simply need to count the 
number of energy levels below the Fermi level 

The eigenvalue problem for the operator /C is similar (although far from being equiv- 
alent in general) to the eigenvalue problem of the particle of unit mass, which is confined 
to move at the positive semi-axis y > and subject to the spike-like potential 

The operator /C is to be compared with the operator U\ of the evolution during the unit 
imaginary time. The latter can be easily diagonalized: 

U x f E = e~ E f E , (A.3) 



with fsiy) = A(y —^f-) , E = E + logA/27r + |g, and A(y) is the modified Airy function 



2 



^ j , ^ ^ , ^ 6 v «/. i 4g , 



A(!j) = I ^-e ipy+ ^ +l £ (A.4) 



where the contour 7 is such that Op 3 > as p — > 00 along 7. The spectrum is determined 

2E 



from the condition that /_b(0) = 0, i.e. A(-W-) = 0. For large values of E this equation 



can be solved using quasi-classics. It gives 



E k - \ \ * - 00. (A.5) 

Assuming that for our problem we may use the same asymptotics we conclude that 

(A0~|^t, (A.6) 

which means that 

a = y. (A.7) 

Remark. It is interesting to note that a similar Schrodinger problem arises in the Born- 
Oppenheimer approximation to the quantum mechanics of a particle in two dimensions 
confined by the potential x 2 y 2 , which is a good model for the matrix potential Tr[X, Y] 2 , 
see [fZIJ. 
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Appendix B. Solution of the one-matrix model from the KP equation and 
double scaling limit. 

As an example illustrating the application of the KP hierarchy to matrix integrals, we 
will derive the critical singularity and the double scaling limit of the one-matrix integral: 



r 3 

Z N (t) = / DM exp iVTr ^ t q M q 

J 0=1 



:b.d 



in the case of a cubic potential. The fact that the matrix integral (|B . 1|) is a tau-function of 
the KP hierarchy has been established in [ 26j . A direct derivation of the Hirota equations 
( |4.24| ) from the matrix integral can be found in P? |. The free energy 



and the specific heat 



u(t) = 2d 2 t F_ 



N 



(B.2) 



(B.3) 



can be obtained from the KP equation ( }4.3ij| ), if we take it into account that the specific 
heat is actually a function of a single parameter: 



u = to 3 ifj 



lit! t 



(B.4) 



and is given in the Gaussian limit by 



t 3 = 0: u= . 



By inserting (|B.4j) into (|4.33|) we derive that the function / obeys the following 
differential equation 



where tf^o is a constant, 



x = — i/j = t$u. 



t| 3£| 



(B.5) 
ordinary 

(B.6) 
(B.7) 



In the large- N limit, we have an algebraic equation for x(i/j): 



dx 



(^-V'o) T - + 2x + 9^ = 0, 
dip 
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whose solution depends on two constants (tpo, f3): 



x 



- — 2 - 3 (V - \i>o 



(B.9) 



By comparing this with the Gaussian limit, we get: i/jq = 0,(3 = — |. Hence the large- N 
result is: 

1 _ 1 

(B.10) 



1 2 2 1 

tits — -to = —ot 3 U 



The critical point is 



■<■<■ =1-1 , tpc 



3u 2 



i 

2 \ 3 

9 



(B.11) 



To compare with the results of Brezin et al. p8[ we set: 



1 9 

tl — 0, to — , t.S — 5" 

' 27V' 6 jyi 



and get x 



— ^-r. Hence 



108^, 



which is exactly the result of the tedious computations |28[ using the distribution of the 
eigenvalues in the large- N limit. Equation ( [4.33| ) also allows a trivial derivation of the 
double scaling limit for the pure 2d quantum gravity [|2U|| . This limit consists in sending 
x—x c to zero and iV to infinity in such a way that the double scaling variable z = N b (x—x c ) 



remains finite. We try the ansatz: ip = ip c + N a v(z). Since ip ~ ip c + const • \fx — x c we 
have: b = —2a. We then plug this ansatz into (|4.33| ) and obtain 



v 



+ 6N 2+5a vv' + -iJ; c N 2+5a + f-zv' + -v) N 2+6a = 0. 
3 V 3 3 J 



:b.i2) 



To keep here the non-linear term (the source of all higher-genus corrections ~ 1/N 2g ) we 
impose the condition: 2 + 5a = 0, which gives 



a 



A 6=*. 

5' 5 



One immediately sees that the last term in ( B.12|) vanishes in the double scaling limit 
(iV 2+6a = iV" — > 0). Integrating ( |B.12[ ) once with respect to z, we finally obtain the 
Painleve I equation: 

v" + 3v 2 = cz, (B.13) 
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where c = —\ib c = — 

J 33 

Remark. We obtained this equation without making use at all of the method of orthogonal 
polynomials - the only technique known until now for these purposes. This raises the 
hope that the method of Hirota equations proposed here will allow us to obtain the non- 
perturbative description (beyond the loop expansion) of some interesting non-critical string 
theories, such as 0(2) model whose grand partition is known to satisfy the KdV 
hierarchy [ p9| . 
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